Dendritic flux avalanches in rectangular superconducting films — numerical simulations 
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Dendritic flux avalanches is a frequently encountered instability in the vortex matter of type II 
superconducting films at low temperatures. Previously, linear stability analysis has shown that such 
avalanches should be nucleated where the flux penetration is deepest. To check this prediction we do 
numerical simulations on a superconducting rectangle. We find that at low substrate temperature 
the first avalanches appear exactly in the middle of the long edges, in agreement with the predictions. 
At higher substrate temperature, where there are no clear predictions from the theory, we find that 
the location of the first avalanche is decided by fluctuations due to the randomly distributed disorder. 

PACS numbers: 74.25.Ha, 68.60.Dv, 74.78.-w 



I. INTRODUCTION 



II. MODEL 



Dendritic flux avalanches have been observed in many 
kinds of type II superconducting films at low tempera- 
tures pp. The origin of the avalanches is a thermomag- 
netic instability mechanism between the Joule heating 
created by vortex motion and the reduction of the criti- 
cal current density as temperature increases [2H1]- 



When transverse applied field is gradually increased, a 
critical state is formed from the edges, with almost con- 
stant current density, and non-zero magnetic flux den- 
sity. The flux penetration is gradual and smooth until 
the conditions for onset of instability is fulfilled. Then, 
an avalanche is nucleated and large amounts of magnetic 
flux rushes in from the edges and forms a complex tree- 
like structure. The onset conditions for such events are 
determined by a competition between the Joule heat- 
ing, heat diffusion and heat removal to the substrate [I]. 
From this theory, it is expected that instabilities are most 
likely to happen where the electric field is highest and the 
flux penetration is deepest. Since rectangles in low ap- 
plied fields experience a flux penetration which is slightly 
deeper in the middle of the long edges [5] , we hence ex- 
pect the first avalanches to be nucleated there. At higher 
fields, the flux front straightens out, and it is less clear 
where the most favorable nucleation location will be. 



In this work, we will run numerical simulations of a 
thermomagnetically unstable superconducting rectangle 
in applied transverse field. The location of the avalanches 
will be discussed in context of prediction from linear sta- 
bility analysis. The simulation method is described in 
Ref. 0. 



Consider a type II superconducting film subjected to 
an applied magnetic field transverse to the film plane, 
ramped at a constant rate H a . The shape of the sample 
is a rectangle with dimensions 2a x 2b. 

A thcrmomagnetic instability is a consequence of the 
nonlinear material characteristics of type II superconduc- 
tors, which conventionally is approximated by a power 
law 
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where E is electric field, J is sheet current, J = |J|, 
po is a resistivity constant, d is sample thickness, J c is 
critical sheet current, and n is the creep exponent. The 
temperature dependencies are taken as 



Jc = Jco(l - T/T c ), n = ni/T, 



(2) 



where T c is the critical temperature. The electrodynam- 
ics must be supplemented by the heat diffusion equation 



cT = kV 2 T- -(T-T q ) 
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where c is specific heat, n is thermal conductivity, h is 
the coefficient for heat removal to the substrate, and To 
is the substrate temperature. 

Eq. ([3| must be solved together with Maxwell's equa- 
tions and the material law, Eq. (IT). The description of 
the simulation method is in Ref. B] and here only the 
key ingredients are outlined. The main challenge is to 
invert the Biot-Savart law in an efficient way. This is 
done by including also the vacuum outside the sample in 
the simulation formalism. At the cost of including the 
extra space, one can use a real space/Fourier space hy- 
brid method with the very attractive performance scaling 
of 0(N log N), where N is the number of grid points. 
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FIG. 1: Distributions of magnetic flux density, B z , at 
H a /J c0 = 0.052(a), 0.064(b), and 0.078(c). Substrate tem- 
perature is To = 0.17T C . The magnetic field is highest at the 
edges, seen as a white rim, while the black central region is 
still flux-free, i.e., B z = 0. The avalanches are all small, most 
of them are fingers, while in (c) some avalanches also have 
two or three branches. 



Quenched disorder is important in order to give re- 
alistic nucleation conditions, and it is introduced as a 
5% reduction of J c0 in randomly selected 5% of the grid 
points. 

The parameters correspond to typical values for MgB 2 
films 0: T c = 39 K, j c0 = 1.2x 10 11 A/m 2 , p n = 7 pflcm, 
k = 170 W/Kmx (T/T c ) 3 , c = 35 kJ/Km 3 x (T/T c ) 3 , 
h = 200 kW/Km 2 x (T/T c ) 3 , and m = 15. Here p n is 
normal resistivity, and we have J c q = dj c0 , and po = p n - 
The sample dimensions are 2a — 8 mm, 2b — 4 mm, 
and d — 0.4 ^m, while the total simulated area is 12 x 
6 mm 2 , where the extra space outside the sample is used 
to implement the boundary conditions. The total area is 
discretized on a 768 x 384 equidistant grid. 




III. ANALYTICAL PREDICTION 

The linear stability analysis in Ref. [5] determines the 
condition for instability onset as 



(4) 



where I* is the threshold flux penetration depth, \J' C \ = 
JcqTo/T c , and E is the background electric field. For a 
given E, the sample is unstable when the flux front ex- 
ceeds I*, which means that avalanches should appear first 
where the flux penetration is deepest. In a rectangle, this 
is in the middle of the two long sides and consequently it 
is expected that the first avalanches appear there. 



IV. SIMULATION RESULTS 

The simulations start from initially zero-field cooled 
state, and the applied field is increased at constant rate, 
PoH a = 10 T/s. The high ramp rate was chosen by per- 
formance reasons, since a high ramp rate closes the gap 
between the velocity of the avalanches and the normal 
flux penetration. Yet, there is clear separation of time 
scales, since full penetration is reached in approximately 
Jco/H a — 6 ms, while the duration of the avalanches is 
less than 0.1 ps. 

Fig. [T] shows the distributions of magnetic flux density 
transverse to the film plane, B z , at H a /J c o = 0.053(a), 
0.064(b), and 0.078(c). The substrate temperature is 
T = 0.17T C , and at this low temperature the thresh- 
old field for avalanche activity is low, H t h — 0.052J c o. 
Image (a) is just above the threshold, where two small 
fingers have appeared symmetrically in the middle of the 
two long edges, just as expected from the theory. In 
image (b) many more avalanches have come at the long 
sides, and for the first time there are avalanches appear- 
ing at the two short sides, also these in the middle. In 
image (c) there are even more avalanches and they now 
cover most of the boundary, except close to the corners. 
All avalanches are small, and will consequently be seen 
as small, jumps towards zero in magnetization curves. 
With increasing field, it seems like a trend of increasing 
avalanche size. 

Fig. [2] shows B z for a simulation run at a higher sub- 
strate temperature, To = 0.22T C . The figure contains 
just one avalanche, which appeared at = 0.12J c0 . 
Both the increased H t h and the much larger size of the 
avalanche is as expected for higher To [E]< Also expected 
is the complex branching pattern [6]. The location is 
not entirely symmetric. The reason is that at deeper flux 
penetrations the flux front straightens out, as seen in reg- 
ular Bean-state of the upper part of the image. Hence, 
the location of the avalanche is instead determined by 
the fluctuations in electric field due to the quenched dis- 
order. The avalanche location is typical for what is seen 
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FIG. 2: Distribution of magnetic flux density, B z , at H a /J c o = 0.14 with substrate temperature To = 0.22T C . Only one 
dendritic avalanche has been nucleated. Note that the location is not in the middle of the side and the avalanche is larger and 
have more branches than those at lower To, in Fig. [I] 

experimentally, e.g., by magneto-optical images in NbN 
rectangles [9]. 

V. CONCLUSIONS 

We have simulated dendritic flux avalanches in super- 
conducting films in the shape of a rectangle. The re- 
sults confirm the prediction from linear stability anal- 



ysis that avalanches will first appear in the middle of 
the long sides, at least for low To. At higher To, the 
avalanche location was not symmetric. In general, the 
results regarding avalanche size, threshold field, and time 
between avalanches follow the typical dependency on T , 
i.e., that increasing T gives larger H t h, larger avalanches, 
more branches, and longer time between the avalanches. 
Hence, the simulations of this work also confirms the cor- 
rectness of the simulation method of Ref . [5] . 
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